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ABSTRACT 

Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for 
angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the labo- 
ratory. A liquid sodium auj dynamo experiment at the New Mexico Institute of Mining and Technology 
provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local 
stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds num- 
ber, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by 
the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability 
conditions are presented and unstable conditions for the sodium experiment are compared with another 
proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number 
achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a 
wide range of wavenumbers and further to observe the transition to the turbulent state. 

Subject headings: accretion, accretion disks — instabilities — MHD — plasmas 



1. INTRODUCTION 

A significant problem in accretion disk theory is the 
nature of anomalous viscosity. In order for accretion to 
occur, angular momentum must be transported outward. 
The central problem in astrophysical accretion disks is that 
observed accretion rates cannot be due to ordinary molec- 
ular viscosity. A robust anomalous angular momentum 
transport mechanism must operate in accretion disks. 

In 1991, the magnetorotational instability (MRI), dis- 
covered by Velikhov (1959) and Chandrasekhar (1960), 
was reintroduced as a mechanism for excitation and sus- 
taining MHD turbulence in a magnetized but Rayleigh- 
stable fluid by Balbus & Hawley (1991a). Since then, many 
numerical and analytic studies of the MRI have been per- 
formed under varying conditions (Balbus & Hawley 1991b; 
Matsumoto & Tajima 1995; Hawley, et al. 1996; Stone et 
al. 1996; Gammie 1996; Sano & Miyama 1999; Noguchi et 
al. 2000; Sano & Inutsuka 2001). Nevertheless, amidst all 
the theoretical attention granted to the MRI, it has never 
been demonstrated in the laboratory. In light of this fact 
an ato dynamo experiment at the New Mexico Institute 
of Mining and Technology provides a unique opportunity 
to study the MRI in a rotating metal annulus using liq- 
uid sodium. In this paper, a local stability analysis is 
performed and the results are compared with theoretical 
analysis from a similar proposed experiment at the Prince- 



ton Plasma Physics Laboratory (Ji et al. 2001; Goodman 
& Ji 2001). Varying aspects of the experiments are dis- 
cussed with stable and unstable regions identified in terms 
of magnetic field strength and shear flow. In addition the 
number of unstable modes and the Prandtl number fur- 
ther define the parameter space. If the number of unsta- 
ble modes is large compared to unity, then there exists the 
possibility of observing turbulence generated by the MRI. 
Finally we investigate the instability boundary when fluid 
turbulence is injected as for example through the Ekman 
layer flow. 

2. NEW MEXICO CtL) DYNAMO EXPERIMENT 

The New Mexico aui dynamo (NMD) experiment is a 
collaboration between the New Mexico Institute of Min- 
ing and Technology and Los Alamos National Laboratory 
(Colgate, et al. 2001). The experiment is designed to cre- 
ate an astrophysical dynamo, the aw-dynamo, in a rapidly 
rotating laboratory system. The apparatus consists of two 
coaxial cylinders (Fig. 1), rotating at different angular ve- 
locities and therefore creating Couette flow in the annular 
volume. Liquid sodium fills the volume between the cylin- 
ders and the end walls. Solid plates attached to and co- 
rotating with the outer cylinder with an angular velocity, 
f^2 define the end walls. (In addition, for the dynamo ex- 
periment, an external source of helicity is supplied, driven 
plumes, but this is not part of the MRI experiment.) The 
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Fig. 1. — A schematic drawing of the New Mexico Institute of Mining and Technology, £7-Phase, of the aui dynamo experiment. The inner 
cylinder is rotated relative to the outer or main cylinder of radius, R2 = 30.5 cm at Q.2 and where the inner cylinder of R\ = 15.2 cm is 
rotated at f2i . The ratio, £7i/f22, can be adjusted from unity to infinity by a brake on f22, not shown, and where Q± is driven. A fixed gear 
ratio also allows the maximum shear case at marginal stability, Hi/Q2 = (R2/R1) 2 = 4 to be obtained. An axial or quadrupole magnetic 
field is supplied by coils shown also schematically. The magnetic field internal to the liquid sodium in the annular space between the two 
cylinders and bounded by the end walls (co-rotating with the outer cylinder) is measured by an aerodynamically shaped probe containing 
six, 3-axis, Hall detectors. These detectors can measure fields from 0.1 G, the earth's field, to up 10 kG. Thus the initial static magnetic field 
distribution, B z , can be measured and compared to the field components, B r and Bg, produced by the MRI. The probe can be mounted 
either radially or axially depending upon the mode information desired. 



schematic of the flow field, (Fig. 1), places particular em- 
phasis on the primary diagnostic of multiple, 3-axis, mag- 
netic field Hall effect detectors (sensitivity: O.f to 10 kG) 
located in aerodynamically shaped probes within the ro- 
tating conducting fluid. We expect that the radial pertur- 
bations from the MRI and their azimuthally sheared result 
will produce a fluctuating B r and Bg field from an orig- 
inal imposed static B z field through MRI growth. These 
fluctuating fields are the result of the linear and non-linear 
growth of the various MRI modes transformed by the dif- 
ference of the sheared Couette flow at a given radius and 
the probe angular velocity, Q 2 , of the outer cylinder. A 
significant difficulty will be the observation of the linear 
growth of any particular MRI mode because the time con- 
stant for establishing the initial axial field within the con- 
ducting liquid sodium will be long, ~ 30/S12, compared 
to the expected growth rate, ~ f^, of the instabilities 
as derived in this paper. We therefore expect to observe 
primarily the near steady state of the non-linear limit of 
various modes, but the sequential linear phases may be 
observed during the comparatively slow rise of the field. 
If the applied field or flux is amplified by the MRI such 
as a dynamo, then we expect to see fluctuating fields sig- 
nificantly greater than the applied field. In addition since 
the inner and outer cylinders are driven separately, the 
relative torque as a function of the applied magnetic field 
becomes an integral diagnostic of the non-linear limits of 
the instability growth. 

By driving the inner cylinder and applying a variable 
brake with a corresponding torque measurement to the 
outer cylinder one can explore the full range of Couette 
velocity profiles including the marginal Couette flow hy- 
drodynamic stability condition discussed next. This con- 
dition of maximum or marginal stable Couette profile can 
be established in the experiment precisely by gear ratios 
and so the degree of turbulence measured by the torque 
can be explored at the stability boundary. In addition the 



pressure will be measured at five radii and compared to 
the pressure distributions expected of the various Couette 
profiles. A finite torque measurement can be interpreted 
in terms of turbulence existing between the two cylinders. 
No turbulence or perfectly laminar flow will exert a torque 
of the order l/R e , R e the fluid Reynolds number where 
Re ~ 10 7 , compared to a turbulent torque, ~ I/i?e 1 / 2 , if 
the Ekman layer circulation leads to the weak turbulence 
that we discuss later. This same possible weak turbu- 
lence can also be measured by introducing a very weak 
field, B min ~ 1 G, small enough so as not to cause the 
growth of MRI in resistive liquid but large enough so that 
an unstable flow or weakly turbulent flow can be measured 
as fluctuations in B r and Bg with the Hall effect probes. 
Therefore the fluid flow conditions can be fully explored 
before the application of magnetic fields designed to create 
the MRI. When the MRI does take place, then the insta- 
bility can be recognized as a departure from the previously 
measured initial fluid state. 

It is critical to have large shear rates in order to observe 
the maximum growth rates of the MRI. However, excessive 
shear will hydrodynamically destabilize the flow by the 
Kelvin-Helmholtz instability. Let us consider a Couette 
flow profile in cylindrical coordinates. Take r, 8, z as the 
radial, azimuthal and axial directions respectively. The 
radial distribution of angular velocity of the flow, fi(r), is 
given by (Landau & Lifshitz 1959) 

_ n 2 R\ - n,Rj 1 (r>! - n 2 )R\Rl 

() ~ R\-R\ r 2 R\-R\ ' [) 
where R\(R 2 ) and f2i (O2) are the inner(outer) radii and 
angular velocities. 

In the limit of infinitely large hydrodynamic Reynolds 
number, R e , the stability condition for Couette flow is 
given by VLiR\ < VL 2 R\ (Landau & Lifshitz 1959). There- 
fore, in order to maximize the shear flow within the appa- 
ratus, the NMD experiment has been designed such that 
RijR\ = 2 and fli/Q 2 = 4, guaranteeing that Sl\R\ = 



3 



VI2 R 2 . I n addition to stability constraints, stress limi- 
tations in the experiment require that an upper limit of 
2 = 33 Hz be placed on the frequency of rotation of the 
outer cylinder. Of course lower rotation rate can be used 
in both NMD experiment and Princeton Plasma Physics 
Laboratory (PPPL) experiment, but it is assumed for this 
analysis that the highest rates are of greatest scientific in- 
terest. 

The assumption of stable Couette flow implies a lami- 
nar flow with no turbulence. On the other hand the initial 
acceleration of the fluid to the final state of Couette flow 
from an alternate initial state implies a transient enhanced 
torque, because, just as in the accretion disk, the laminar 
friction is too small. However, in the experimental ap- 
paratus, the transient, Couette flow profile is Helmholtz 
unstable so that turbulence is a natural and expected re- 
sult of the " spin- up" of the flow. 

On the other hand the Ekman flow creates a relative 
torque between Oi and VL2 that we expect to be bal- 
anced by a weak turbulence as observed by (Taylor 1936) 
and analogous to the spin-up turbulence. This turbulence 
may also influence the stability conditions but primarily 
the ability to distinguish turbulence caused by the MRI 
from the hydrodynamic turbulence caused by the Ekman 
layer. We therefore analyze the MRI stability conditions 
as a function of hydrodynamic turbulence preexisting in 
the liquid and therefore of the Prandtl number. At large 
enough levels of turbulence, the effective electrical resistiv- 
ity can also be increased and therefore decrease the mag- 
netic Reynolds number, R m , and therefore influence the 
conditions of excitation of the MRI. 

In comparison to the New Mexico Dynamo Experiment, 
a similar experiment at the Princeton Plasma Physics Lab- 
oratory has been proposed to look for the MRI in a rotat- 
ing liquid metal annulus (Ji et al. 2001). The PPPL exper- 
iment utilizes liquid gallium, an easy to handle metal with 
properties similar to liquid sodium (see Table 1). Note, 
however, the higher density and higher resistivity of liquid 
gallium, which limit the maximum rotation speed and the 
maximum achievable R m . The dimensions of the PPPL 
experiment are slightly different, enabling them to acquire 
larger shear flow rates, Ji et al. (2001). For Ri = 5 cm 
and i?2 = 15 cm then R2/R1 = 3 with a typical Vti/Q 2 = 
9. The conditions for instability for both experiments are 
discussed in Sec. 4. 

3. LOCAL STABILITY ANALYSIS 

The angular velocity of Couette flow confined between 
coaxial cylinders with radii Ri < r < R2 and cylindrical 
angular velocities fii, is given by 

n(r)=a+A (2) 

where we define a and b as 

O2-R2 — Vl\R 2 



a = 



b = 



(fil — Q2)R\R\ 



(3) 



The incompressible and dissipative MHD equations de- 
scribing the dynamics of liquid metals are given as fol- 
lows, 

V • B = 0, (4a) 



V - V = 0, 
dB 



dt 



= V x (V x B) +r]V A B, 



dV ,„ ^ \ t 7" (B-W)B 1 / B 2 



(4b) 
(4c) 



(4d) 



where B is the magnetic field, V is the velocity, 77 is the 
magnetic diffusivity, p is pressure and v is the kinematic 
viscosity. In cylindrical symmetry the system of equations 
(4a)-(4d) have stationary solution Vq = (0,rQ(r),0) and 
Bq = (0, Bgoif), B z q), where B z q is a constant, Bgo oc 1/r, 
and the angular velocity profile, fi(r), is given by expres- 
sion (2). One needs to investigate the time evolution of 
perturbations to this equilibrium state governed by the 
linearization of system of equations (4a)-(4d). A sim- 
ilar analysis of perturbations was performed by Good- 
man & Ji (2001), who showed that local WKB approxi- 
mation gives results for growth rates of instability, which 
are close to the growth rates obtained by the solving full 
boundary value problem in radial direction. As it is es- 
pecially stressed by Ji et al. (2001) and Goodman & Ji 
(2001), WKB local analysis leads to a good approxima- 
tion to the growth rates even in the case of the scale of 
perturbations being equal to or comparable to the sizes 
of the vessel. Thus, in this work, we limit ourselves to 
the local approach, which is much easier to carry out than 
the full eigenmode analysis, because it allows one to ob- 
tain an algebraic dispersion relation. The perturbations 
b = (b r ,bg,b z ), v = (v r , v$,v z ) are assumed to be axisym- 
metric and proportional to exp(7i — ik z z — ik r r) where 7 
is the associated growth rate. It is also assumed that the 
minimum possible wave numbers in r and z directions are 
Kmin = tt/ (i?2 — Ri) and k zm i n = tt/L. The linearized 
equations of motion are then given by 



0=1 ik r j 17 //,--_ 



0=1 ik r I b r — ik z b z , 



jb r — —ik z B z0 v r — rjk 2 b r , 
dfl , 

769 = -tk z B z0 vg + — — b r 



-TV, 



d\nr 
d_ (Boq 
dr 



T]k bg, 



7« r — 2Qvg = — 



- — ik z B z0 b r H b t 

Airp \ r 

p 



iv z = 



4^o y kzBz ° be ~ br 

-vk 2 vg, 

ik z B z0l , ., pi 



d_ 

dr 



b z +ik z — -vk 2 v z , 



where the epicyclic frequency k is defined as 
2 1 d(r 4 n 2 ) .„„ dn 2 



dr 



AVI 2 



dlnr ' 



(5a) 

(5b) 
(5c) 

(5d) 
(5e) 

) 

(5f) 
(5g) 

(6) 
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Table 1 

Actual and Normalized Quantities in the Sodium and Gallium Experiments 



Property 


Actual 
Sodium Gallium 


Normalized 
Sodium Gallium 


Kinematic Viscosity, v[cm s ) 


7.M0~ 3 


3.2T0~ 3 


3.6-10" 8 


2.2-10-'' 


Reynolds Number, R e 


- 


- 


1.3 • 10 7 


3.0- 10 6 


Magnetic Diffusivity, r](cm 2 s _1 ) 


810 


2000 


4.2T0" 3 


1.4T0" 1 


Magnetic Reynolds Number, R m 






120 


4.7 


Density, p(g cm -3 ) 


0.92 


6.0 






Alfven Speed, VA(cm s _1 ) (10 3 gauss) 


2.9 • 10 2 


1.1 • 10 2 


4.6-10~ 2 


1.2T0" 1 


Inner Radius, i?i(cm) 


15.25 


5 


.5 


.33 


Outer Radius, i?2(cm) 


30.5 


15 


1 


1 


Length, L(cm) 


30.5 


10. 


1 


0.66 


Inner Angular Velocity, Sli (s _1 ) 


829 


533 


4 


8.2 


Outer Angular Velocity, f^s -1 ) 


207 


65 


1 


1 


Prandtl Number, P M = R m /R e 






9.2 • 10~ 6 


1.6 • 10~ 6 


Ekman Turbulent Prandtl =ff, Put 






0.012 


6.3 x 10~ 3 



pi is the perturbation of the pressure, and k = \/k 2 + k 2 
is total wave number, respectively. Note that k can be 
expressed through a and b in Eq. (3) as 



k 2 = 4a I a + 



(7) 



which vanishes when velocity shear is maximum (a = 0, 
b = r 2 n). 

These equations lead to the following local dispersion 
relation 

r 2i 2 k 2 

(7 + vk 2 )(~i + rjk 2 ) + k 2 V A 



U2 



+K ^ 1 + ri k 2 y + -— k/v A 



rflnr 



+ 

where 



2ik z V ASz 



0, (8) 



Va6z = 



BanB 



eo n zo 



B 



Va = TZZ (9) 



4irp ' ' 71 4irp 

Neglecting all 1/r terms compared to k yields the follow- 
ing dispersion relation, which is identical to the dispersion 
relation derived by Ji ct al. (2001) 



[{ 1 + vk 2 )( 1 + r 1 k 2 ) + (k z V A ) 2 ] 2 

+ K 2 {j + T]k 2 ) 2 + 



rlO 2 

-'■ 1 (k z V A ) 2 



k 2 

= 0. (10) 
= 0, and hence, 



dmr 

In the case of maximum shear flow, a 
k = 0, the dispersion relation simplifies to 

[( 7 + ^ 2 )(7 + V k 2 ) + (k z V A ) 2 ] 2 - 4fi 2 ^ = 0, (11) 
which immediately yields the following solutions for 7 



7= 2 



-{v + ri)k 2 



±\ (v + v) 2 k 4 - 4 vqk A + k 2 z V 2 ± 



2rik 2 V A 
k 



,(12) 



Only when we take the plus sign for the square root term 
and the minus sign for the last term in eq. (12), does it 
give the unstable solution, and all the other three solutions 
are stable. 

The MRI occurs only when the second term inside the 
square root of (12) is negative, i.e., 



LJ^LU V +UJ A <2 



k z fluj A 



(13) 



where 



vk 2 



rjk 2 and uj a = V A k z . Thus, vis- 



cosity, magnetic diffusion and magnetic tension stabilize 
the MRI, whereas the shear flow destabilizes it. The con- 
dition for neglecting v and 77 can be derived from (12) by 
evaluating the expression under the square root, i.e., 



[y - ry) 2 fc 4 « 



4uj A (2k z Q — ku A ) 
k ' 



For k r — 0, eq. (14) further reduces to 



{ri-vfkl ^ 4(2n-u; A ) 



LO A 



(14) 



(15) 



Thus, it is apparent that magnetic diffusivity and kine- 
matic viscosity only affect high k z modes. 

If the magnetic diffusivity is large and the applied mag- 
netic field is weak, eq. (12) reduces to two roots 7 = —vk 2 
and 7 = —rjk 2 . The former root corresponds to a hydro- 
dynamical branch in which the fluid is disconnected from 
the electromagnetic force and behaves as a pure fluid. The 
latter root is an electromagnetic branch, in which the mag- 
netic field diffuses as in vacuum. With an increasing mag- 
netic field, bifurcations occur in which the hydrodynamic 
and electromagnetic branches are split into four branches. 
When uj v ~ lo a <C uj v the unstable solution of (11) is given 

by 

LU A (2fl — U> A ) 



7 



-oj v + 



(16) 



showing that the unstable solution emerges from the hy- 
drodynamical branch. Though magnetic diffusion dimin- 
ishes the MRI, the branch remains unstable if the condi- 
tion lo a > uj v oj u k/ (2k z Q) is satisfied. Notice however, that 
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even if the magnetic diffusivity is high, a weak magnetic 
field is capable of generating the MM. 

Next, let us consider the hydrodynamical limit(V^ = 0) 
with arbitrary Couette flow profiles, (a ^ 0). In this case, 
k =/= in general and we have to go back to eq. (10) for 
deriving solutions. Two trivial solutions are 7 = —-qk 2 , 
corresponding to the electromagnetic branch and the sec- 
ond is given by 

-y = -uk 2 ±i K ^, (17) 
k 

which is the hydrodynamical branch. 

The effect of finite n and high v is shown in Fig. 2, where 
the rotation speed of the cylinders, viscosity, and magnetic 
diffusivity of Fig. 2a corresponds to the point C of Ji et 
al. (2001), and the wavenumber is fixed at (k z , ry) = (1, 1) 
for Fig. 2a, 2c and (4, 1) for Fig. 2b, 2d, respectively. The 
growth rates of the four roots of eq. (10) are shown as 
a function of the axial magnetic field strength B z . The 
epicyclic frequency n is finite in Figs. 2a and 2c, whereas 
k = in Figs. 2b and 2d. The kinematic viscosity v is 
taken as the actual value of Gallium [Fig. 2a] and Sodium 
[Fig 2b], whereas we make it artificially high {v > 77) in 
Figs. 2c and 2d to see the effect of anomalous increase of 
v due to possible turbulence. 

In the hydrodynamical limit (B z = 0) in Fig. 2a, the so- 
lutions of the hydrodynamical branch, (7 ~ 0) are complex 
[eq. (17)]. Near B z = 2000 gauss, these solutions are sep- 
arated and both become real. Only one solution becomes 
unstable. Thus, if the flow is not a maximum shear flow 
profile (k 0), the MM is stabilized for weak magnetic 
fields. 

Figures (2c) and (2d) show that the turbulence sup- 
presses the unstable MM mode. Since the Ekman layer 
may make the fluid weakly turbulent, it is important to es- 
timate the scale of this turbulence, which will be discussed 
in Sec. 5. 

In the next section, stability diagrams are presented and 
compared for both experiments. 

4. STABILITY DIAGRAMS AND GROWTH RATES IN 
SODIUM AND GALLIUM EXPERIMENTS 

The comparison between the NMD and Princeton ex- 
periments is done by comparing their typical parameters 
in Table 1. In order to evaluate the physical differences 
between the experiments, we compare the dimcnsionless 
parameters, presented in the second column of Table 1. 
We use R 2 and fi^ 1 as units of length and time to obtain 
the dimensionless quantities. 

We choose fii = 84.8 Hz and Q 2 = 10.34 Hz as the 
typical values for the gallium experiment, which corre- 
sponds to point C of Ji et al. (2001). The global magnetic 
Reynolds number becomes 

Rm = R 2 Q 2 (R 2 -R 1 ) 

V 

The R m is higher in the sodium experiment, which is de- 
signed to observe the au dynamo (Colgate, et al. 2001; 
Pariev 2001). All the growth rates are evaluated at the 
radius r = f, which satisfies O(f) = V^i^2- We also 
define the global fluid Reynolds number as 

-Z?2^2(i?2 — i?i) 



R e 



(19) 



Using these parameters, the growth rate is obtained 
by solving eq. (10) numerically, and the unstable regions 
are plotted in Fig. 3 as a function of axial magnetic 
field strength and wave number k z . The unit of k z is 
n/(R 2 - Ri), k r is ir/L in Figs. 3, 4, 5, 6, 8, 9 and 10. The 
minimum possible values for the dimensionless k r and k z 
are unity. We fixed k r as unity in Fig. 3. Notice that in 
both experiments, a strong field suppresses high k modes 
because of the magnetic diffusivity and magnetic tension. 
In the sodium case (Fig. 3a), higher k z modes are desta- 
bilized, and the growth rate is higher compared to the 
gallium case (Fig. 3b). In the gallium case, the suppres- 
sion of the unstable modes with low magnetic field occurs 
due to finite k (see Sec. 3 for details). 

Figures 4 and 5 demonstrate the dependence of the 
growth rate on the wave numbers in the sodium and gal- 
lium experiments. In Fig. 4, an axial magnetic field B z is 
fixed at 3 x 10 3 gauss, and in Fig. 5, at 400 gauss. When 
B z = 3 x 10 3 gauss, a number of k z modes(fc z < 8) are 
destabilized in the sodium experiment (Fig. 4a), while in 
the gallium experiment (Fig. 4b) only the k z = 1 mode is 
destabilized. 

A more significant difference between the sodium and 
gallium experiments is shown in Fig. 5. For weak magnetic 
fields no mode is unstable in the gallium experiment (Fig. 
5b), whereas higher k z modes (k z > 50) are excited in the 
sodium experiment (Fig. 5a). As we noted before, the finite 
k suppresses the unstable MRI modes with weak magnetic 
field. 

In order to show the importance of the maximum shear, 
we plot the growth rate for the gallium experiment with 
k = in Fig. 6. We reproduce Figs. 3b, 4b and 5b in 
Figs. 6a, 6b and 6c respectively, except we take 2 = Oi/9 
for the maximum shear [see Eqs. (3) and (7)]. In Fig. 6a, 
high k z modes are unstable with low magnetic field, which 
are stable with finite ft(Fig. 

rcffig3b). Maximum velocity shear also destabilizes high 
k r modes (Figs. 6b and 6c), so many modes will be excited 
in the gallium experiment at maximum shear. In all cases, 
the {k Zl k r ) = (1, 1) mode is dominant, but mode coupling 
may occur in the nonlinear regime to excite turbulence in 
the gallium experiment. 

While the maximum shear flow profile leads to easy ex- 
citation of MRI, it falls on the border line for pure hy- 
drodynamical instability. We show the contour plot of 
the unstable region for the modes {k r ,k z ) — (1,2) and 
(3, 5) of the sodium experiment in Fig. 7. In the regime 
Q. 2 /Q,\ < 0.25, sodium is hydrodynamically unstable. The 
maximum shear flow is indicated as a solid line. High 
wavenumber modes are unstable only near the maximum 
shear flow and weak magnetic field in the hydrodynami- 
cally stable region (Fig. 7b) . Growth rates for the finite k 
shear case are shown in Fig 8. High wave number modes 
are stabilized compared to the k = case, Fig. 6, but sev- 
eral modes are still unstable. Comparison of Fig. 3b to 
Fig. 8a and Fig. 4b to Fig. 8b shows that even for highly 
sub-critical flow with fl\/fl 2 = 2 the sodium experiment 
will allow one to observe more MRI unstable modes than 
the gallium experiment for slightly sub-critical flow with 
fii/fia = 8.2 (point C of Ji et al., 2001). Thus, the sodium 
experiment has a higher potential for the observation of 
turbulence due to the nonlinear development of the MRI 
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Fig. 2. — The complex stability plane of the combinations of limiting Couette and non-limiting Couette flow profiles with laminar and 
turbulent flows. Figs, (a) and (c) are for a non-limiting Couette flow profile (the gallium experiment), and Figs, (b) and (d) are for the 
limiting Couette flow profile (the sodium experiments). In Figs, (a) and (b), the flow is assumed as laminar(f <C rj), and (c) and (d), the 
flow is turbulent, v > r\. The corresponding Prandtl numbers are Pm* = 1-20 for (c), and 3.23 for (d). A remarkable difference between the 
non-limiting Couettc(a) and limiting Couette (b) profiles is the hydrodynamical branch, where {—vk ) is stable with the non-limiting Couette 
profile and weak magnetic field, whereas it becomes unstable with a weak magnetic field and the limiting Couette flow profile. With strong 
turbulence, [Figs, (c) and (d)], the non-limiting Couette flow profile is always stable (c), whereas the limiting Couette flow profile becomes 
unstable with a strong magnetic field (B z > 2 X 10 3 gauss). In Figs, (a) and (c) (k z , k r ) = (1, 1); in Figs, (b) and (d) {k z , k r ) = (4, 1). 
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Fig. 3. — Growth rates of the MRI in sodium (a) and gallium (b) experiments Growth rates are shown as a function of fc z and B z . All 
the contours of contour figures are equally spaced, and only positive growth rate contours are shown. The radial wavenumber k r is fixed to 
be unity for both calculations. When B z = 3 X 10 3 gauss, fc z = 4 mode is the most unstable mode with 7 = 300s -1 , and for the sodium 
experiment (a), whereas k z = 1 with 7 = 25s -1 for the gallium(b). Higher k z modes arc destabilized with weak magnetic field in the sodium 
experiment, but with strong magnetic field, they are suppressed by the magnetic diffusivity and magnetic tension. Weak magnetic field modes 
in gallium are suppressed because of finite n (see Fig. 6) 



at the rotation profiles with smaller shear, i.e. when the 
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(a) 
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kz 

Fig. 4. — Growth rates as a function of wavenumbers. An axial magnetic field strength is fixed at 3000 gauss for both experiments. The 
value of the maximum contour is 350s -1 for (a) and 100s _1 for (b). The growth rate decreases with k r for both experiments. Up to k z = 8 
modes arc unstable in the sodium case (a), whereas only the k z = 1 mode is unstable in the gallium case (b). Other parameters are taken 
from Table 1. 
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Fig. 5. — Same as Fig. 4 with B z = 400 gauss. The value of the maximum contour is 130s 1 for (a) and 7.0s 1 for (b). No mode is 
destabilized in the gallium case (b), whereas high k z modes are excited in the sodium case (a). 



flow is highly stable in the absence of the magnetic fields. 

5. TURBULENCE DERIVED FROM THE EKMAN LAYER 

FLOW 

We believe we need to understand the minimum ex- 
pected turbulence level in the fluid before attempting to 
observationally separate the MM growth from an un- 
known background. Of course that background will be 
measured first, but this section deals with an estimate of 
its magnitude. First we must point out that if there were 
no initial turbulence in the flow, then the expectation of 
laminar flow is that the torque or power required to drive 
the experiment would be negligibly small. Instead, since 
we know that there must be a torque associated with the 
Ekman layer flow and therefore a minimum power required 
to drive the Couette flow, we calculate this and compare 
it to the observed instability based upon the gradient of 
angular momentum (Richard & Zahn 1999). 

Richard & Zahn (1999) have extensively reviewed the 
earlier experimental work on Couette flow of Wendt (1933) 



& Taylor (1936). In the maximally stable case where the 
outer cylinder rotates at and the inner one is stationary, 
VL\ = 0, and thus where the angular momentum increases 
outward, the flow is observed to be weakly unstable to a 
finite amplitude instability despite the prediction of sta- 
bility for a positive angular momentum gradient. This 
observed instability has been invoked by Richard & Zahn 
(1999) as a possible mechanism for the a-viscosity of Ke- 
plerian accretion disks. The implication is that the free 
energy of the flow is accessed through a finite amplitude 
instability producing turbulence. However, this turbulence 
in the positive angular momentum gradient, theoretically 
stable regime, is observed to be much weaker then the 
inverse case, exerting a much smaller torque than the tur- 
bulence generated by the unstable Couette flow, which 
in turn is self-excited by the flow of free energy through 
the turbulence itself. This lack of equivalence causes us 
to ascribe the occurrence of turbulence in the stable flow 
case to the back-reaction to the torque of the Ekman layer 
flow. We calculate this back reaction turbulence here as 
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Fig. 6. — Growth rates of the MRI in the gallium experiment with re = 0. All the parameters are the same as the point C in Ji et al. 
(2001) except Q.2 = Cli/9. The value of the maximum contour is 35s -1 for (a), 100s —1 for (b) and 30s _1 for (c). (a): Growth rate with 
k r = 1. When B z < 4.5 X 10 3 gauss, k z = 2 and higher modes are destabilized due to maximum shear rate (re = 0). (b): With B z = 3.0 X 10 3 
gauss. High k r modes, which are stable in the finite re case (Fig. 4b), arc unstable, (c): With B z = 4.0 X 10 gauss. High wavenumber modes 
(k r < 15, k z < 22) are unstable, whereas no mode is unstable in the finite re case (Fig. 5b). 




Fig. 7. — The MRI unstable region of sodium experiment for two different modes. Mode number is fixed as (k r , k z ) = (1, 2) (a) and (3, 5) 
(b). The value of the maximum contour is 270s -1 for (a) and 250s -1 for (b). The region ^2/^1 < 0.25 is hydrodynamically unstable. 
Parameters of point A correspond to those of Figs. 3a and 4a, B to those of Fig. 5a, C to Figs. 8a and 8b, D to Fig. 8c, respectively. 



the likely lower limit of turbulence for these experiments. 
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Fig. 8. — The MRI growth rates with finite k flow shear, Q.2/Q.1 = 0.5, for sodium experiment (Table 1). (a) The wavenumber k r is fixed 
as unity; (b)With B z = 3 ■ 10 3 gauss; (c)With B z = 1.5 X 10 3 . The value of the maximum contour is 170s for (a) and (b) and 65s _1 for 
(c). High wavenumber modes are suppressed in each case. 



The initial experimental measurement of the torque as well 
as the in situ magnetic fluctuations for very low fields can 
and will be compared to these predictions. Finally we 
note that the predicted ratio between the torques for the 
stable and unstable flows, based on the Ekman layer flow 

1 /2 

G s tabie/G uns tabie — Re /R e , roughly agrees with the mea- 
surements of Taylor (1936). We delay until a later paper 
a full analysis of this problem, but give here an estimate 
of this expected turbulent viscosity as compared to the 
purely laminar one. 

An Ekman layer forms adjacent to the surfaces of the 
end walls. This flow is both radial and azimuthal, thin, 
laminar, and high speed. The resulting flux of angular 
momentum creates a torque on the fluid and an unstable 
velocity profile between the inner and outer cylinders. The 
unstable shear flow at both the inner and outer cylindrical 
boundaries results in a "law of the walls" or "logarith- 
mic profile" turbulent boundary layer (Schlichting 1960) 
with each of the cylinder walls. This turbulence extends 
from the inner to the outer differentially rotating cylinders 
and creates the turbulent stress necessary to transfer the 
torque between them. 



The assumption of the stability of Couette flow in the 
experiment is limited by the formation of an Ekman layer 
adjacent to the surfaces of the end plates. Since these end 
plates corotate with the outer frequency ^2, then at any 
radius r < R 2 the fluid will be rotating faster than the end 
wall. An Ekman layer forms (Prandtl 1952), when the cen- 
trifugal force is not balanced by a pressure gradient. The 
pressure in the Ekman layer is the same as the pressure in 
the bulk of the cylinder but the centrifugal force is smaller 
in the Ekman layer because of friction with the end wall. 
As a result, a (negative) radial flow develops in a thin layer 
of thickness S with a mean radial velocity < v r >~ rO/2 
while undergoing a mean azimuthal motion < vg >~ rO/2 
. The analysis of Prandtl (1952) results in the thickness 
6 ~ rj\fW e = J v/0,2 — 5.6 x 1CP 3 cm in sodium and 
7.0 x 10 -3 cm in gallium, using the parameters of Table 1. 
Hence a radial (negative) current, F r , flows of order 

F r ~ -(5(27ri? 2 )i?2^2/2 = -ttR^/^/rI em's" 1 . (20) 
This (negative) radial flow at both ends towards the 
axis must be balanced by a positive, slower radial flow 
throughout the central region. The Ekman flow merges 
with the central flow by a boundary layer at the inner 
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cylinder surface resulting in a circulation within the Cou- 
ette flow volume driven by the Ekman layers at each end. 
Since this flow represents a flux of angular momentum, 
pF r i?2^2/2, from the inner radius, Ri, to the outer ra- 
dius, R2, there must be a torque, G, transmitted by the 
fluid corresponding to the difference in the flux of angular 
momentum between these two surfaces or 

G E = pF r {R\ - i??)0 2 /2 = p(3ir/8)Rlnl/\/R e . (21) 
where we have used the ratio R1/R2 = 1/2 for the sodium 
experiment. (A factor of 5.9 x 10~ 3 smaller is implied for 
the gallium experiment.) Not only does this torque, be- 
tween Ri and R2 determine the power required to drive 
the flow, but also imposes a requirement for a weak tur- 
bulence within the so-called stable Couette flow in order 
to transmit this torque between i?i and i?2- This level 
of turbulence becomes the minimum effective viscosity or 
turbulent viscosity, v tl of the MR1 experiments. By way 
of comparison in addition we calculate the torque as if 
the flow were completely laminar and compare the two 
torques. 

The shear stress for turbulent or laminar flow, r t , tl, 
is characterized by either an effective turbulent viscosity, 
v t or laminar, v^. In the turbulent case the angular mo- 
mentum flux from the Ekman layer must be balanced by 
the viscous stress from the rate of shearing, A = rdQ/dr = 
—2VL2R\t 2 [see Pringle (1981)] resulting in a viscous drag 
per unit area, Tt = pvtA = — 2v t £l2R\r~ 2 p and therefore 
a torque per unit length, t t = —2nr 2 T t = —A-KVt^R^P- 
This has the proper scaling since the torque must be in- 
dependent of radius. The corresponding laminar torque 
per unit length, = —2irr 2 TL = —Attvl^R^P where the 
viscosity, ul, is fixed by the fluid properties and not vari- 
able with the strength of the turbulence. Then the total 
laminar torque per half length becomes, Gl = = 
2Tny L n 2 LR 2 p = 2TrLRjnf/R e . The ratio of the two 

torques becomes Gl/Ge = (16/3)(L/i?2)-Re • Since 
R e is very large and L/R 2 = 1, the laminar torque is 
negligibly small compared to the Ekman layer torque and 
therefore torque balance requires turbulence to enhance 
the effective viscosity. This effective turbulent viscosity is 
obtained by equating the Ekman angular momentum flux 
to the viscous shearing torque giving 

v t = (3/16)^^£^i? 2 = 10.0 & 12.7 cm 2 a' 1 , (22) 

for the sodium and gallium experiments respectively. 

The structure of this turbulence is problematic. It has 
been described by (Taylor 1936), as initially a series of long 
parallel vortices, "Taylor columns" that extend the full 
length of the annular space and that at greater R e these 
columns break up becoming fully developed turbulence. 
We expect at some value of magnetic field strength that 
these vortices will be suppressed by magnetic field of suf- 
ficient strength. It seems unlikely, however, that a return 
to laminar flow would take place, because the flow profile, 
with no turbulent shear stress, but still the Ekman flow, 
will become more distorted from the stable profile result- 
ing in stronger turbulent drive. However, for now we defer 
analysis and expect guidance from future experiment. 

Finally a turbulent magnetic Prandtl number, Pm«, for 
measuring the strength of this Ekman driven turbulence, 
can also be estimated as 

P Mt = v ± = 0.012 & 6.3 x 10~ 3 (23) 
V 



respectively. Despite the very small size of the Ekman 
layer, the turbulence generated by such a flow influences 
the ability to distinguish turbulence caused by the MRI at 
low values of magnetic field from the hydrodynamic tur- 
bulence caused by the Ekman layer. At higher values of 
the magnetic field, above that affected by the Ekman tur- 
bulence, the effects caused by the MRI should be clearly 
recognizable. 

Figure 9 shows the dependence of the MRI growth rate 
with the turbulent v t as a function of an axial magnetic 
field strength and the Prandtl number for the most un- 
stable modes in both experiments. In sodium experiment, 
higher k r modes are also unstable. Figure 10 shows the 
growth rate for (k r , k z ) = (3, 4) mode. Note that this mode 
is stable in gallium, even with laminar viscosity only. Fi- 
nite k for the gallium experiment prevents the MRI from 
developing with weak magnetic field(compare Figs. 3b and 
6a, Figs. 4b and 6b above), and MRI exists only in the re- 
gion 3 x 10 3 gauss < B z < 6x 10 3 gauss (Fig. 3b). However, 
the sodium experiment with maximum shear, k = 0, can 
be destabilized with already very weak magnetic field, in 
the range of 50 gauss < B z < 100 gauss. 

In conclusion, the presence of Ekman layers is signifi- 
cant for the determining of the power necessary to sustain 
the differential rotation in the apparatus but has a negli- 
gible effect on the condition of the excitation of the MRI 
as it has been already mentioned briefly in Ji et al. (2001). 
However, the turbulence excited due to the presence of Ek- 
man layers may interfere with our measurements of per- 
turbations of magnetic field excited by the MRI. There- 
fore, the presence of weak turbulent perturbations due to 
Ekman layers seems unavoidable whenever one observes 
the excitation of the MRI. It also seems unlikely that the 
imposed magnetic field in both sodium and gallium exper- 
iments can significantly exceed a value of a few thousands 
Gauss. The characteristic amplitude of the perturbations 
of the magnetic field due to the Ekman layer turbulence is 
~ B z Xv t /ri ?s PutRz- The typical value of such perturbed 
magnetic fields is of order of 1 % of the applied field (see 
Table 1 for Put), thus, limiting the possible MRI mea- 
surements of growing fields to more than 1 % of the initial 
magnetic field. 



6. DISCUSSION AND CONCLUSIONS 

There are several aspects of each experiment that war- 
rant discussion. The first deals with the analysis per- 
formed for both experiments. In this paper we used the 
full dispersion relation which depended not only on an az- 
imuthal magnetic field but also has all terms proportional 
to 1/r retained. This corresponds to the geometrical effect 
of the curvature from the cylindrical geometry. Only when 
one neglects all 1/r terms in the dispersion relation does 
one obtain the results of Ji et al. (2001). This allows us 
to consider many different magnetic field configurations, 
some of which will be suitable for studying the MRI. 

Nevertheless, even with these significant problems and 
differences, it should be noted that both the NMD and 
PPPL experiments have an excellent chance of observing 
the MRI in the laboratory. Both experiments obtain very 
high growth rates under varying conditions yielding a flex- 
ible set of opportunities. 
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Fig. 9. — The effect of turbulence on the MRI for the sodium (a) and gallium (b) cases. The growth rate is plotted as a function of an axial 
magnetic field and a magnetic Prandtl number. Wavenumbers arc fixed for the maximum growth rate of each experiment [(k r ,k z ) = (1,4) 
for sodium, (1, 1) for gallium]. The value of the maximum contour is 280s -1 for (a) and 23s -1 for (b). The growth rate decreases when the 
turbulent viscosity increases. The MRI can be sustained in fully turbulent fluid in both cases for the minimum Ekman layer driven turbulence, 
(P m = 0.012 for sodium, P m = 6.3 X 10~ 3 for gallium). 
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Fig. 10. — The effect of the turbulence on the MRI for the mode (k r ,k z ) = (3,4) in the sodium experiment. The value of the maximum 
contour is 90s -1 . Note that this mode is stable in gallium even without turbulence. This mode will remain unstable even with turbulence. 
The critical Prandtl number for this mode to become stable is P m ~ 0.3. 



Finally we note that the effect of an azimuthal magnetic 
field and an analysis of nonaxisymmetric modes are still 
open problems. Noguchi ct al. (2000) showed that the local 
dispersion analysis in shear flow may fail even for the qual- 
itative estimation of growth rates. The cigenmode analysis 
for nonaxisymmetric modes is necessary for further under- 
standing of the MRI instability in the NMD experiment. 
We are now developing the shooting method code for solv- 
ing Eqs. (4a)- (4d) simultaneously. Spatial dependence of 
the radial wavenumber and azimuthal magnetic field de- 
pendence will be analyzed. 
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